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Abstract 

An alternative perturbative expansion in quantum mechanics which allows a full expression of 
the scaling arbitrariness is introduced. This expansion is examined in the case of the anharmonic 
oscillator and is conveniently resummed using a method which consists in introducing an energy 
cut-off that is carefully removed as the order of the expansion is increased. We illustrate this 
technique numerically by computing the asymptotic behavior of the ground state energy of the 
anharmonic oscillator for large couplings, and show how the exploitation of the scaling arbitrariness 
substantially improves the convergence of this perturbative expansion. 
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1. Introduction 

Perturbation theory, as everybody knows, consists of calculating the quantities associated to 
a Hamiltonian H(g) as Taylor expansions in the coupling g around g = 0, which corresponds to an 
exactly solvable Hamiltonian. In systems like the anharmonic oscillator, however, these series fail 
to converge. This problem is generally addressed by resumming the series with the Pade or Borel 
methods. In cases like the double well potential or physically interesting field theories like QCD, 
however, these methods are inapplicable due to the existance of instantons. 

In this paper, we alternatively look at perturbation theory in a different, very physical way. 
The central idea is the realization that all the energy scales of the unperturbed Hamiltonian partic- 
ipate in the calculation of an observable (an eigenstate, for instance) of the perturbed system at a 
certain energy. This view naturally leads us to investigate two aspects of perturbation theory. On 
one hand we investigate how every energy scale contributes to the perturbed eigenstate and how 
the calculation of these contributions can be optimized using energy-dependent scaling transforma- 
tions. This leads us to introduce a perturbative expansion that exploits the scaling arbitrariness to 
a maximum. On the other hand, we identify the divergences in the anharmonic oscillator as due to 
the fact that the contributions calculated in perturbation theory from arbitrarily high energies are 
unbounded. From this realization we propose a resummation technique based on the introduction 
of an energy cut-off. 

The optimization of perturbation theory using scaling properties is accomplished in Quantum 
Field Theory (QFT) by the Renormalization Group (RG) 1 ' 2 . In the process of renormalization 
a new mass scale fj, needs to be introduced and we get vertex functions expressed in terms of 
more parameters than the physically independent ones. This parameter redundancy is exploited 
by the RG: we can suitably change the scale jj, and the rest of parameters in such a way that 
we always remain in the same physical theory. Although these changes do not affect the theory, 
they do modify its perturbation expansions order by order. Thus, depending on the energy scale 
at which we are working, we can select the appropriate fj, scale which optimizes this perturbative 
expansion. As is well known, this technique has proved to be extremely fruitful in QFT, its most 
important achievement being, perhaps, the discovery of asymptotic freedom in non Abelian Gauge 
Theories 3 ' 4 . 

However, it is not clear to us that the conventional use of the RG fully exploits the scaling 
properties of the theory. In the computation of vertex functions at a certain energy, we will 
obviously use a \i scale of that same energy. In that computation, however, all scales are involved 
since we perform loop integrals over the full range of virtual momenta; the effect of the perturbation 
on these other scales should, nevertheless, be optimally computed using a different fi. 

In a recent paper 5 , we showed how the RG in a 4-dimensional (ft 4 theory could be extended 
to provide improved perturbation expansions free of infrared divergences in the massless case. 
There, we computed higher order vertex functions (n > 4) using skeleton expansions in which full 
propagators and 4-point functions are inserted respectively at each line and vertex of the skeleton 
diagrams. In this way, depending on the momenta running through a line or reaching a vertex of a 
diagram, the propagators and couplings conveniently adjust themselves to those particular scales. 
This would take care of the objection of the previous paragraph for higher order vertex functions. 
Unfortunately, the equivalent treatment for the primitive divergences would imply solving the 
Schwinger-Dyson integral equations 6 , which is a highly non-trivial matter. 

In this paper we provide, for the simple case of one-dimensional Field Theory (i.e., Quantum 
Mechanics), a perturbative expansion that maximally exploits the arbitrariness of scaling in the 
sense stated above. That is, we describe a perturbative expansion in which, in the computation at 
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a certain energy, the effect of the interaction on all the other energies involved is suitably computed 
using the optimal scale for every energy. 

Consider the Hamiltonian H(g) = P 2 /2m + V(g,x) where g = corresponds to the exactly 
solvable system and g = gf to the perturbed one. We shall consider the eigenvectors and eigenvalues 
of H{g) as functions of g. However, instead of doing the usual Taylor expansion in g, we will use 
another expansion which is much more convenient for visualizing how the various energy scales 
contribute to the calculation and for exploiting the scaling arbitrariness. We will divide the interval 
g G [0,5/] into N equal segments by introducing N — 1 intermediate points and will transport all 
eigenvectors and eigenvalues from one point to the next using first order perturbation theory. N is 
here the order of our expansion; when N is increased, the intermediate segments become smaller 
and smaller and first order of perturbation more and more accurate. 

In this framework, scaling can be easily implemented to improve perturbation theory at every 
step for each energy level. If gi and g 2 are the initial and final points of a certain elementary step, 
for every given energy level we appropriately rescale H(g 2 ) (U~ 1 (X)H(g2)U(X) = X 2 H(g2), where 
U(X) is the scaling operator given in coordinate representation by U(X)^(x) = \/~Xty(Xx)) in such 
a way that H{gi) and H(g 2 ) are as similar as possible at that particular energy. This requirement 
can be expressed by demanding that the classical turning points of the two potentials coincide at 
that energy, since this ensures that the corresponding wave functions will be essentially confined to 
the same space interval, minimizing their difference. We then compute the eigenvector \I> of -ff (52) 
using first order perturbation theory and finally determine the eigenvector \E' of H(g2) applying 
the scaling operator to ^ = U{X)i>). 

We are essentially dividing the interaction introduced at each step into two parts: an inter- 
action that does not change the natural scale of the eigenvectors and which is computed in first 
order of perturbation theory, and a pure change of scale which is computed exactly by applying the 
appropriate scaling transformation. In this way, we remove from perturbation theory the burden of 
having to express the changes of scale introduced by the interaction, and it can thus be more accu- 
rate. It is important to recognize that these changes of scale introduced by the interaction strongly 
depend on the energy and, therefore, the scaling transformation U{X) will vary from eigenvector 
to eigenvector. 

Having introduced our alternative expansion and its optimization by scaling in section 2, we 
turn in section 3 to the discussion of the divergences of perturbation theory, particularizing for the 
anharmonic oscillator. 

As is well known, the divergence of the Taylor expansions for the energy levels of the anhar- 
monic oscillator is attributed to the fact that for negative couplings the potential is unbounded from 
below and the system becomes unstable. This instability manifests itself in the form of singularities 
in the second and third sheets of the Riemann surface associated to the analytic continuation of 
the eigenvalues for complex couplings 7,8 . These singularities have g = as accumulation point and 
hence, the perturbation expansions for the energy levels have a null radius of convergence. The 
knowledge of the analytic properties of the eigenvalues 7,8 and of the asymptotic behavior of the 
Taylor coefficients 9 enables us to use the powerful machinery of the Borel Transformation 10,11 to 
resum the series. 

In the framework of the perturbation expansion presented in this paper, we look at the di- 
vergence in a different way. For large N, one would expect that the errors made by first order of 
perturbation theory in an elementary step are of order (Ag) 2 = (gj/N) 2 ; thus, the total error for 
the N steps would go like 1/N making the series converge to the right result. However, this is not 
true for the anharmonic oscillator. In this case, the perturbation changes the asymptotic behavior 
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of the potential for large x and, consequently, the perturbation becomes arbitrarily large for high 
energy states (which are sensitive to the potential at large x). Thus, the errors made by first order 
of perturbation theory in an elementary step cannot be uniformly bounded with respect to energy 
by increasing N. The divergence in the low energy states then arises because these big errors at 
high energy are propagated down to the lower energy states in the subsequent elementary steps. 

The use of the optimizing rescalings explained above ameliorates this problem reducing sub- 
stantially the errors at high energy. Nevertheless, although the difference between the perturbed 
and unperturbed potentials effective at each energy is greatly reduced by making their classical 
turning points coincide, this difference is still arbitrarily large at high energies and the divergence 
persists. 

We show how this problem can be overcome introducing an energy cut-off by, for example, 
limiting the number of states of the system. In this way, we can reduce the errors of first order 
of perturbation theory for all states as much as we want by increasing the number of elementary 
steps. The theory can then be resummed by carefully removing this cut-off as we increase the 
order of perturbation. 

Thus, we regard our technique as a double expansion: in the number of elementary steps (the 
order) and in the number of eigenstates that we take around the energy at which we are computing 
(the cut-off). The way in which we calculate the double limit of this expansion is crucial: it can 
lead to the correct result or to a divergence. We provide a criterion to find the suitable order 
of perturbation associated to every value of the cut-off and show how, as expected, this order is 
substantially lower for the expansion that fully exploits the scaling arbitrariness. 



2. An alternative perturbative expansion and scaling 

In this section we will introduce the mentioned alternative expansion and will show how to fully 
exploit the scaling arbitrariness to optimize it. We have in mind its application to the anharmonic 
oscillator, but the formalism can be presented in a general way: we will consider a Hamiltonian of 
the form H(g) = P 2 /2m + V(g,X) where the potential satisfies V(g, —x) = V(g,x), V(g,x) < 
V(g,x') if \x\ < \x'\ and lim^ x \^ (yo V(g,x) = co. 

If ^ n (g), E n (g) are respectively the eigenvectors and eigenvalues of H (g), we will first deter- 
mine the expressions of their derivatives with respect to g in terms of them. The defining equations 
are: 



The eigenvectors are completely determined up to a global phase factor. This phase can 
be conveniently chosen to satisfy the equation (ty n (g), 4f n (g)) = (3). Differentiating (1) with 
respect to g yields: 



H{g)^ n (g) = E n {g)^ n {g) 
(* n (g),* n (g)) = 1 



(1) 
(2) 



H(g)*M + H(g)* n (g) = E n (g)* n (g) + E n (g)* n (g) 



(4) 
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And using (3) and the projection of (4) on ^k(g) Vfc, we get te desired expressions: 



E n (g) = (* n (g),H(g)* n (g)) (5) 

h(g),H(g)V n {g 
E n (g)-E k (g) 



By repeatedly differentiating these equations with respect to g and setting g = 0, we obtain the 
ordinary perturbation theory. As we pointed out in the introduction, we will use here a different 
strategy which is much more convenient for later introducing the scaling arbitrariness. We shall 
alternatively expand in the number of intermediate points in which we divide the interval g £ [0, g/} 
(the 'order of perturbation') and eigenvectors and eigenvalues will be propagated from one point 
to the next using a first order approximation. In view of (5)-(6), this first order approximation 
should be: 



^ E n (g)-E k (g) 
E n {g + A<?) ~ £ n (s) + (*„(<?), (tf ( 5 + Ag) - #(<?))<M<?)) (8) 

but it is, nevertheless, more convenient to replace (7)-(8) by 

¥„(» + Ag) = (*„( 9 ) + J«„(s))/||*„(s) + Jtf„( 9 )|| (7') 

where ».(,) = E (» t te).(gto + ^)-Hte))».te)) fftto) 
^ E n {g)-E k {g) 

E n {g + A<?) ~ (* n ( 5 + A 5 ), + Ag)* n (g + A 5 )) (8') 

because, although these formulas coincide with (7)-(8) at first order, they have the advantage 
of yielding exactly normalized states and the exact energies corresponding to the (approximate) 
eigenstates for finite values of Ag. 

No scaling properties have yet been introduced; so far we have just changed the way in which 
we do our expansions. The different philosophies underlying each of these two expansions, can 
probably be best illustrated in the more familiar context of expanding a function around a point. 

Suppose that we want to solve the equation 

^(x,y) = O v g(x,y) (9) 

with the boundary condition g(0,y) = f(y), (O y is an operator that acts on the y variable only). 
Here are the two strategies: 

i) Differentiating (9) with respect to x and setting x = 0, we get: 



0(O,») = O^(O,») = O»/(y) 
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and we can build the Taylor expansion 

N 



N—>oo ' re! 

n=0 

ii) We divide the interval [0, x] into N parts and transport g(£, y) from one point to the next using 
the first order approximation 

9{i + e, y) ~ <?(£, 2/) + eO y s(£, y) (e = x/iV) 
By repeatedly doing this, we get: 

g(e,y)~f(y) + eO y f(y) 
5 (26, y) ~ y(e, y) + eO y y(e, y) = /(y) + 2eO y f(y) + e^/ (y) 

5 (7Ve, y) ~ g((N - 1)6, y) + eO y g((N - l)e, y) = £ ^ e n O™/ (y) 

n=0 ^ ' 

Thus, using this approach we get the expansion 

g(x,y)= lim V- — r O"/(y) (11) 

yv jv-oo-^ (N-n)\N n n! y v ; v 7 

n=0 v 7 

To apply this to the expansion of a function around a point, we just need to realize that the 
function g(x,y) = f(x + y) is determined by the equations 

|^(x,y) = g(x,y) ; y(0,y) = /(y) 

which is the former situation with O y = d/dy. Substituting x by x — x$ and y by x in (10) and 
(11) we get the expressions 



t) f{x) = lim Yl " ( 12 ) 

N^oo £ — ' re! 

JL ATI f^txn) 

ii) f{x) = lim Y — Ar J - { -^l( x -x ) n (13) 

Term by term, (12) and (13) coincide as N — ► oo, but not uniformly so, and as a result (13) 
has certain advantages over (12); for example, for the function f(x) = 1/(1 + x), (12) converges 
in the interval [0,1] while (13) converges in [0, a] where a is the solution of log(a) — 1/a = 1, 
(a ~ 3.59). We therefore may wonder whether our expansion based on (7'), (8') already has 
improved convergence properties over conventional perturbation theory. 

However, let us go further to improve our expansion by the exploitation of the scaling ar- 
bitrariness. As is well known, the eigenfunction corresponding to a given energy level, quickly 
vanishes outside the classically allowed region. Thus, the distance between the classical turning 
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points of the potential is a measure of the length scale characteristic of that eigenstate. This scale 
can be preserved for all couplings by doing a scaling transformation which depends on g. 

The scaling operator is given in the coordinate representation by 

U(\)V(x) = Vx^(Xx) (14) 
and transforms the position and momentum operators in this way: 

U- 1 (X)XU(X) = X/X U- 1 (X)PU(X) = XP (15) 
The Hamiltonian transformed by the operators U(X(g)) is then 

H(g) = X- 2 (g)U-\X(g))H(g)U(X(g)) = ^ + X~\g)V (g , X / X(g)) (16) 

and X(g) can then be chosen to maintain the scale corresponding to the nth eigenstate by solving 
the equations 

E n (0) = V(0,x ) E n (0) = X- 2 (g)V(g,x /X(g)) (17) 

The first equation determines the classical turning point xo of the unperturbed potential for 
the energy E n (0) and the second equation says that X(g) should be chosen by demanding that 
the classical turning point for the perturbed potential at the same energy E n (Q) is again xq. The 
coupling g can, of course, be conveniently reparameterized to make sure that equal increments of 
the coupling correspond to similar increases of the interaction in the rescaled Hamiltonian H. 

The perturbative expansion for the nth state associated to H(g) will clearly be better than 
that associated to H{g): ^f n (gf) and ^ n (0) have the same number of nodes and are confined to 
the same region. In this way (^ n (0), *& n (gf)) is maximized and the terms (^(O), ^> n (gf)) (k ^ n), 
which express the strength of the perturbation, are minimized. The eigenvector and eigenvalue 
corresponding to H(gf) can then be immediately obtained from ty n (gf), E n (gf) by 



* n (g f ) = U (X(g f))^> n (g f) (18) 
E n (g f ) = X 2 (g f )E n (g f ) (19) 

as follows from (16). 

What it has been done so far is just setting the appropriate scale for the nth eigenstate. This 
would correspond in Field Theory to set the \i scale equal to the energy at which we are calculating. 
At this stage, both the ordinary Taylor expansion and our alternative expansion could be applied 
to the rescaled Hamiltonian H(g). From here on, we will suppose that we are studying a particular 
eigenstate and that we have already done the mentioned global rescaling of the Hamiltonian and 
dropped the tildes from the rescaled quantities. 

As one can see in (7) or (7'), all energy scales participate in the calculation of a certain 
eigenstate and the effect of the perturbation on each of these energies, should be calculated at their 
appropriate scales. It is in taking care of this fact where our alternative perturbative expansion 
is much more convenient. Its advantage is that it merely consists in repeating many times a very 
simple routine: first order of perturbation theory. Thus, the scaling arbitrariness can be expressed 
in a single step very simply, and this new algorithm repeated for every step. Instead, in a Taylor 
expansion, the iVth order is something very difficult to interpret. 
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Obviously, the way to exploit the scaling arbitrariness at a single step is just to do, for every 
eigenstate and for the elementary increment of the coupling, the global rescaling that we carried 
out for our selected state and the entire range of the coupling. 

To be more precise, suppose that gi and g 2 are the initial and final points of a certain ele- 
mentary step and that we know ^i(gi), Ei(gi) V7. The appropriate rescaling of H(g 2 ) for the fcth 
state is 

H{g 2 ) = X- 2 U- 1 (X k )H(g 2 )U(X k ) (20) 
where, in analogy with (17), X k is determined by 

E k (gi) = V( gi ,x k ) E k ( 9l ) = XfV(g 2 ,x k /X k ) (21) 

We now use a first order approximation between H(g 1 ) and H{g 2 ) followed by a scaling 
transformation between H(g 2 ) and H(g 2 ). Using (7) and (8) we have the approximations 

*k(92)**k(9i) + S* k (9i) ( 22 ) 

E k (g 2 ) * E k ( 9l ) + (* k ( 9l ), (H(g 2 ) - H( gi ))* k ( gi )) (23) 
which will then be inserted in 

^(52) = U(\ k )V k (g 2 ) (24) 

E k (g 2 ) = \ 2 k E k (g 2 ) (25) 

In the spirit of the more convenient formulas (7') and (8') we have 

*fc(52) ^ (ttfc(si) + 5y k ( 9l ))/\\y k ( 9l ) + 5* k ( gi )\\ (26) 
which we then substitute in 



^>k(g2) = U{\ k )^ k (g 2 ) (24) 
E k (g 2 ) = ^ k (g 2 ),H(g 2 )^ k (g 2 )) (27) 

Since U(X k ) is unitary, ^ k {g 2 ) is exactly normalized. 

Thus, for each state, the approximations are applied to an interaction that maintains the scale 
at that particular energy and the change of scale is carried out by an exact scaling transforma- 
tion. We have, therefore, accomplished our goal: in this perturbative expansion, the effect of the 
interaction at every energy is calculated at its appropriate scale. 
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3. The anharmonic oscillator: cut-off resummation and scaling. 



In this section, we will investigate the improvements in the perturbation theory of the an- 
harmonic oscillator when the scaling arbitrariness is exploited as described in the former section. 
An additional difficulty is, however, present in this case: the instability of the system for negative 
couplings gives rise to a singularity at g = 0. We will show that this singularity is not overcome 
by the exploitation of scaling and that can be interpreted as due to the fact that for any order, 
the perturbative approximations are inadequate at very high energies. To overcome this problem, 
we will introduce a cut-off by limiting the number of states considered and will provide a criterion 
for finding the appropriate order of perturbation that should be associated to every value of the 
cut-off to conveniently resum the series. We will finally show how this order is substantially lower 
in the case in which the scaling arbitrariness is fully exploited. 

We will choose the ground state to carry out all the numerical analysis. An analogous devel- 
opment could be done for any eigenstate. We first perform the global rescaling appropriate for the 
ground state explained in section 2. 

H(g)=P 2 + X 2 +gX 4 ([X,P]=i) (28) 



H(g) = \- 2 (g)U- 1 (X(g))H(g)U(X(g)) (16) 

Using equation (17), we can immediately find that the relationship between a and [3 is 

a = l- [3 (29) 
Then, the rescaled Hamiltonian (which is naturally parameterized with (3) is 

H(J3) = P 2 + (1- (3)X 2 + PX 4 (30) 
/3(g) and \(g) are the solutions of 

(3 = g/\ 6 /? 3 + (5 -2_3 )/3 2 + 3/3 _ 1 = 

1-/? = 1/A 4 A=(l-/3)" 1 / 4 
and E (g) is finally determined by 

E (g) = X 2 (g)E ((3(g)) (32) 

where Eq(/3) is the eigenvalue of the rescaled Hamiltonian (30), which will be calculated pertur- 
batively. It should be noted that while g ranges from zero to infinity, 0(g) goes from zero to one. 
This fact illustrates very clearly the importance of this rescaling of the Hamiltonian in calculations 
for large couplings. From (31), we have that, for large g 



0(- 



g 2/3 3g 4/3 ^ 



\*(g)=g 1/3 (l + 



3g 2 / 3 



(33) 
(34) 



By inserting these equations in (32) we can determine the asymptotic behavior of E (g): 



em = ,'/»{&(!) + - ^(D)^r + o(^r)} (35) 

We will later on calculate the coefficient E (l). This is obviously the most unfavorable quantity 
to calculate in a perturbative expansion since it is the most removed from the non interacting theory. 

We will study our alternative perturbative expansion for the Hamiltonian (30) in the two 
versions presented in section 2. Expansion I will be the one that does not take advantage of 
scaling; eigenvectors and eigenvalues are propagated from one point to the next using (7)-(8) or 
(7')-(8'). Expansion II will be the one that fully exploits the scaling arbitrariness, propagating each 
eigenvector and eigenvalue using the (step and energy dependent) optimal rescalings described in 
(22)-(27). 

As we pointed out at the beginning of this section, we will need to limit the number of states of 
the system in order to resum the series. By limiting the number of states, we mean that we consider 
the projection of the Hamiltonian on the space generated by a finite number of the eigenstates of 
the harmonic oscillator. Since the Hamiltonian is invariant under parity, the parts corresponding 
to the even and odd states are uncoupled. Thus, only even states will need to be considered in the 
calculation of the ground state energy. 

Expansion I can be readily applied to this truncated Hamiltonian. In fact, it is easy to see that 
only a finite number of states contribute to a given order of perturbation and, therefore, beyond 
this number, the truncated expansion gives the exact result at that order. For example, only 8 
states contribute at 3rd order and 26 at 4th order. 

The limitation of the number of states, brings about a slight inconvenience in expansion II 
because the scaling operator does not map our finite dimensional subspaces into themselves. To 
overcome this difficulty it should be noted that everything after equation (20) in the last section, 
formally holds if we substitute U(Xk) by an arbitrary unitary operator. Thus, the obvious thing to 
do is to replace U(Xk) by unitary operators that leave invariant our finite dimensional spaces and 
coincide with U (Afc) when the number of states is taken to infinity. Such operators can be obtained 
by taking the projection of f7(A&) on the finite dimensional spaces in the basis formed by our finite 
number of (even) eigenstates of the harmonic oscillator and unitarizing it by orthonormalizing 
Gramm-Schmidt the columns of this projection starting from the one corresponding to the lowest 
eigenstate up. 

For a large number of states (large cut-off), these operators will virtually coincide with U(Xk) 
at low energy. There will be, however, some differences for the highest energy states which are 
necessary to accommodate the cut-off. Fortunately, as it will be shown, this boundary effect 
introduced by the cut-off does not spoil the scaling optimization. 

The two different ways of carrying out the first order approximations in both expansions I and 
II that were presented in section 2, give rise to two versions of each expansion. We will denote by 
A those expansions in which the states are not exactly normalized at each step (eqs. (7)-(8) and 
(22)-(25))and by B those in which they are (eqs. (7')-(8') and (26), (24), (27)). 

The even eigenstates of the harmonic oscillator will be denoted by | I). In coordinate repre- 
sentation are given by 

ll) = ( ^(2l)\ ) 1/2H2l{x)e ~ X2/2 l = ^- ((P 2 + X 2 )\l) = (Al + l)\l)) (36) 
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We will say that the cut-off is equal to n when we limit the number of states to the set 
I 0),| n). 

The numerical calculations that follow were done using 'Mathematica' on a NeXT computer. 
The programs involved are quite simple: one just needs to write the algorithm for transporting 
eigenvectors and eigenvalues from one intermediate point to the next, and ask the machine to 
repeat this algorithm N times. 

In tables I and II, we list the values of Eo(g = 1) computed respectively with the A versions 
of expansions I and II as functions of the order of perturbation and the cut-off. The exact value is 
1.39235 12 . 



Table I (expansion 1(A)) 



I cut-off/ order— ► 


1 


2 


3 


4 


5 


6 


2 


1.4671 


1.4177 


1.4060 


1.4025 


1.4009 


1.4000 


4 


1.4671 


1.4177 


1.4102 


1.4003 


1.3957 


1.3950 


6 


1.4671 


1.4177 


1.4482 


796.31 


3.0 x 10 9 


9.3 x 10 2U 


8 


1.4671 


1.4177 


1.4556 


3.3 x 10 4 


7.7 x 10 10 


5.7 x 10 4e 


10 


1.4671 


1.4177 


1.4556 


4.0 x 10 e 


3.5 x W 2b 


1.0 x 10 es 



Table II (expansion 11(A)) 



I cut-off/ order— > 


1 


2 


3 


4 


5 


6 


2 


1.4671 


1.4178 


1.4061 


1.4025 


1.4009 


1.3999 


4 


1.4671 


1.4175 


1.4035 


1.3994 


1.3976 


1.3965 


6 


1.4671 


1.4175 


1.4038 


1.3995 


1.3976 


1.3965 


8 


1.4671 


1.4175 


1.4039 


1.4087 


1.3982 


1.3965 


10 


1.4671 j 


1.4175 


1.4038 


1.8093 


4.9 x 10 s 


3.5 x 10 ia 


12 


1.4671 


1.4175 


1.4037 


7.3734 


1.2 x 10 iU 


-5 x 10 27 



It is clear that when we remove the cut-off by taking it to infinity, we get divergent series 
in both cases. Thus, the exploitation of scaling in expansion II is not enough to overcome the 
divergence although, as one can see by comparing tables I and II, it tames it considerably. 

To better understand this divergence and the effects of exploiting scaling, we have listed in 
tables III and IV (for expansions I and II respectively) the energies and norms of all states at every 
intermediate point for the case cut — off = 6, order = 5, g = 1 (/? = 0.43015 . . .). 
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Table III (expansion 1(A)) 



points^ 





1 


2 


3 


4 


5 


energy(O) 


1 


1.021 


1.034 


1.043 


7.710 


2.2 x 10 y 


norm(O) 


1 


1.0005 


1.0008 


1.183 


1.1 x 10 4 


5.4 x W L2 


energy (1) 


5 


5.623 


5.991 


6.448 


2224 


7.3 x 10 11 


norm(l) 


1 


1.026 


1.047 


11.71 


2.0 x 10 5 


3.4 x 10 i3 


energy (2) 


9 


11.25 


13.27 


27.73 


5346 


1.9 x 10 11 


norm (2) 


1 


1.187 


1.844 


24.59 


1.0 x 10 5 


1.3 x 10 ia 


energy (3) 


13 


17.92 


29.19 


573.7 


8.8 x 10° 


7.9 x 10 ia 


norm (3) 


1 


1.629 


7.237 


716.3 


2.2 x 10 6 


1.7 x 10 14 


energy(4) 


17 


25.62 


75.17 


7013 


4.6 x 10° 


3.6 x 10 i4 


norm (4) 


1 


2.395 


20.52 


557.1 


4.7 x 10 B 


1.0 x 10 i4 


energy (5) 


21 


34.35 


166.4 


1791 


4.1 x 10 5 


5.7 x W L2 


norm (5) 


1 


3.434 


10.65 


152.9 


5.6 x 10 5 


5.5 x 10 ia 


energy (6) 


25 


44.12 


234.8 


5697 


1.0 x 10 y 


3.3 x 10 i3 


norm (6) 


1 


2.911 


18.36 


812.2 


1.6 x 10 a 


1.2 x 10 14 



Table IV (expansion 11(A)) 



points^ 





1 


2 


3 


4 


5 


energy(O) 


1 


1.021 


1.034 


1.043 


1.050 


1.055 


norm(O) 


1 


1.0005 


1.0008 


1.001 


1.001 


1.001 


energy (1) 


5 


5.501 


5.843 


6.115 


6.344 


6.542 


norm(l) 


1 


1.002 


1.002 


1.002 


1.003 


1.003 


energy(2) 


9 


10.61 


11.63 


12.42 


13.08 


13.65 


norm (2) 


1 


1.004 


1.005 


1.005 


1.005 


1.005 


energy(3) 


13 


16.21 


18.11 


19.67 


21.12 


22.56 


norm (3) 


1 


1.007 


1.009 


1.012 


1.015 


1.019 


energy(4) 


17 


22.38 


26.80 


30.27 


34.00 


37.96 


norm (4) 


1 


1.019 


1.062 


1.076 


1.091 


1.112 


energy(5) 


21 


32.41 


43.00 


51.39 


60.50 


71.54 


norm (5) 


1 


1.217 


1.418 


1.458 


1.495 


1.530 


energy(6) 


25 


51.09 


98.46 


166.1 


233.3 


295.5 


norm(6) 


1 


1.300 


1.521 


1.541 


1.538 


1.528 



Let us look at the first point. The amount in which the norm of the states differ from 1 
is a measure of the errors made by perturbation theory in the first step. These errors increase 
with energy and are substantially smaller for expansion II showing very clearly the remarkable 
improvement introduced by scaling. 

As one can see in table III, the divergence in expansion I occurs because the big errors at high 
energy are amplified and propagated to the lower energy states in the subsequent steps, finally 
reaching the ground state. The same thing happens in expansion II but for a larger value of the 
cut-off (c± 10). The errors of perturbation theory at a given energy are reduced by increasing the 
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order but then, we have more steps and the errors from higher energy states better propagate down 
to the ground state. Thus, the obvious thing to do is to keep the cut-off fixed while we increase 
the order to reduce all errors as much as we want. 

Table V lists the values of E (g = 1) calculated with expansion 1(A) for a cut-off equal to 6 
as we increase the order of perturbation. 



Table V (expansion 1(A)) 



order 


E (g = i) 


order 


E (9 = 1) 


order 


E (g = l) 


1 


1.4671 


7 


1.0 x 10 32 


13 


1.3934 


2 


1.4177 


8 


4.7 x 10 38 


14 


1.3933 


3 


1.4482 


9 


5.2 x 10 36 


15 


1.3933 


4 


796.3 


10 


1.3937 


20 


1.3930 


5 


3.0 x 10 9 


11 


1.3936 


25 


1.3929 


6 


9.3 x 10 20 


12 


1.3935 


30 


1.3928 



At low orders, as we already pointed out, only the lowest energy states give substantial con- 
tributions to E and we get decent approximations. From orders 4 to 9 there are enough steps 
to let the big errors at high energy propagate down to the ground state: the numbers take off. 
Finally, beyond 10, the order is high enough to curb the errors at all energies originating excellent 
approximations. 

As we pointed out in the introduction, we are considering a double expansion: in the cut-off 
and in the order; and the way in which we take these quantities to infinity is crucial. If we first 
remove the cut-off we get a divergence, while by taking the order to infinity first and then the 
cut-off, such divergences are avoided. 

Table III shows how the norms of the intermediate states grow very rapidly causing the energies 
(which are given by eq. (8)) to take off as well. Thus, the exact normalization of the states in the 
B versions of the expansions should represent a significant improvement. Table VI lists the values 
of E (g = 1) at order 6 as we increase the cut-off, calculated with expansion 1(B) and seems to 
indicate that there is no divergence in this case. 



Table VI (expansion 1(B)) 



cut-off E (g=l) 


cut-off E (g = l) 


1 1.39819 

2 1.39654 

3 1.39351 

4 1.39299 

5 1.39312 

6 1.39338 


7 1.39323 

8 1.39333 

9 1.39339 

10 1.39337 

11 1.39340 

12 1.39339 



Whether the B series without cut-off are in fact convergent or not is somewhat irrelevant 
because, as one can see in table VI, the best approximation is obtained for a finite cut-off (4 in 
this case). This is due to the fact that for a fixed order, the approximations of perturbation theory 
at each step, even though greatly improved and curbed by normalization, are not valid for high 
energy states. 
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As we said, the resummation procedure that we are proposing is to change the order in which 
we take the limits in our double expansion: first take the order to infinity and next the cut-off. 
However, for a fixed cut-off, it is only worth increasing the order until the precision obtained is of 
the order of the error made by the fact that we are limiting the number of states. Next, we will 
give a criterion to determine when we should stop increasing the order for a fixed cut-off. 

In order to do that, we should first know how fast the series for fixed cut-offs tend to their 
limits when we increase the order N. In the A expansions, both the states and energies are 
propagated at each step with first order of perturbation theory; thus, the error at each step 
is of order (Ag) 2 oc 1/N 2 and the total error will go like 1/N. The same applies, in the B 
expansions, to the eigenstates, but not to the ground state energy. This is calculated with the 
formula E = (ty,H*S>). The (approximate) ground state \? can be written as Vl/ = ^q + x(1/N)/N, 
where is the exact ground state and x(l/iV) has a finite limit when N — > oo. Thus, 

E = (tt ,#*o) + ^{(*o,H X (0)) + (x(0),iTtto)} + 0(l/N 2 ) (37) 

Since minimizes Eq, the term proportional to 1/N in (37) must vanish and we get that the 
errors in the B expansions go like 1/N 2 for large N. This fact shows once again the superiority of 
the B expansions. 

If Eq(N) — E (oo) goes like 1/N for the A expansions and like 1/N 2 for the B ones, then 
Qo(N) = (E {N) - E (2N))/(E (2N) - E (AN)) must respectively tend to 2 and 4 for large N. 
In tables VII and VIII we numerically verify this behavior in the particular case cut-off=4, g = 1. 
Something completely analogous is found for the A and B versions of expansion II since the use of 
scaling does not modify the arguments presented. 



Table VII (expansion 1(A)) 



order(N) 


E (N) Qo(N) 


50 


1.392672024 1.98 


100 


1.392527424 1.99 


200 


1.392454389 


400 


1.392417667 


Table VIII (expansion 1(B)) 


order(N) 


E (N) Q (N) 


50 


1.392387789 4.10 


100 


1.392382513 4.04 


200 


1.392381227 


400 


1.392380909 



We can now give a criterion for when we should stop increasing the order for a given cut-off. 
As we said, this should be that the precision obtained with that order is about the error made by 
the fact that we are limiting the number of states. 
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Let Eo(n, N, g) be the ATth order approximation to E${g) when the cut-off is set equal to n. 
From the knowledge of the asymptotic behavior we can estimate the precision obtained. For the 
B versions of the expansions it will be 

e = E (n, N, g) - E (n, oo , g) ~ ^ (E (n, N, g) - E (n, 2N, g)) (38) 

On the other hand, a good estimate of the error made by the fact that we have a finite cut-off 
is the difference A = E (n — 1, N, g) — E (n, N,g) (39). Thus, a sensible criterion would be, for 
example e < A/2, i.e. 

Eo(n-l,N,g)-Eo(n,N,g) > |(£b(n, N, g) - E (n,2N, g)) (40) 
and the appropriate approximation for a cut-off equal to n would finally be: 

E (n,oo,g) = E {n,N,g) - e ~ ±(4E (n,2N,g) - E (n,N,g)) (41) 

where ./V is an order that satisfies (40). 

As an example, we will apply this technique to compute the coefficient E (l) of equation (35) 
that determines the asymptotic behavior of the ground state energy for large g and which, as we 
mentioned, is the most unfavorable quantity of the anharmonic oscillator to calculate perturba- 
tively. In table IX we list the smallest order N for which (40) is satisfied and the corresponding 
approximation E (n, 1) = (4E (n, 2N, 1) — E Q (n, N, l))/3 for a cut-off ranging from 3 to 6 and for 
both expansions I and II (B versions). 



Table IX 





expansion 1(B) 


expansion 11(B) 


n 


N 


E (n,l) 


N 


E (n,l) 


3 


6 


1.0650 


6 


1.0648 


4 


7 


1.0614 


8 


1.0612 


5 


19 


1.06046 


13 


1.06044 


6 


100 


1.060389 


42 


1.060389 



The coefficient Eq(1) has been computed by variational methods 8 , the result being 1.060362 . . . 
It is remarkable the high degree of accuracy that is obtained for a relatively low cut-off. For 
instance, for n = 6 the error made is about 3 x 10~ 5 . This is because we have taken advantage of 
scaling: the free and interacting ground states have the same scale and are, therefore, very similar; 
the small difference between them can then be very well approximated by a linear combination of 
a few higher energy harmonic states. 

Finally, it should be noted how, as expected, expansion II (which makes full use of scaling) is 
much more efficient than expansion I. This is more notable the higher the cut-off is, since it is at 
high energies where the optimizing rescalings are larger. As an example, we compare in Table X 
the values of E (8, N, 1) for expansions I and II. For instance, we only need to calculate to order 
40 in expansion II to get the same precision as with order 100 in expansion I. 



15 



Table X 





expansion 1(B) 


expansion 11(B) 


N 


Eo(8,N,l) 


E (8,N,1) 


20 


1.0611724 


1.0605232 


40 


1.0605931 


1.0604097 


60 


1.0604749 


1.0603898 


80 


1.0604311 


1.0603830 


100 


1.0604105 


1.0603799 



Summary 

We have introduced an alternative perturbative expansion which is very convenient for inves- 
tigating the contributions of the various energy scales in the calculation of an eigenvalue and which 
can be very naturally improved by carrying out energy-dependent rescalings. We have identified 
the divergences in the expansions for the anharmonic oscillator as due to the fact that the errors 
made by perturbation theory in this case cannot be uniformly bounded for all energies. To over- 
come this problem, we have introduced an energy cut-off giving rise to a double expansion in both 
the order and the range of energies considered. Finally we have provided a criterion for how this 
limit should be taken in order to get the correct physical results. 
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